STAR FORMATION IN THE LOCAL UNIVERSE 

Editors Corinne Charbonnel and Thierry Montmerle 

Proceedings of the Evry Schatzman School held in Aussois, September 2010 



STAR FORMATION IN MOLECULAR CLOUDS 

Ralf S. Klessen 1 

Abstract. Stars and star clusters form by gravoturbulent fragmenta- 
tion of interstellar gas clouds. The supersonic turbulence ubiquitously 
observed in Galactic molecular gas generates strong density fluctua- 
tions with gravity taking over in the densest and most massive regions. 
Collapse sets in to build up stars. 

Turbulence plays a dual role. On global scales it provides support, 
while at the same time it can promote local collapse. Stellar birth is 
thus intimately linked to the dynamical behavior of parental gas cloud, 
which governs when and where protostars form, and how they contract 
and grow in mass via accretion from the surrounding cloud material. 
The thermodynamic behavior of the star forming gas plays a crucial 
part in this process and influences the stellar mass function as well as 
the dynamic properties of the nascent stellar cluster. 

This lecture provides a critical review of our current understand- 
ing of stellar birth and compares observational data with competing 
theoretical models. 



1 Overview 

When wc look at the sky on a clear night, we can note dark patches of obscuration 
along the band of the Milky Way. These are clouds of dust and gas that block the 
light from stars further away. Since about one century we know that these clouds 
are the birthplaces of new stars. With the current set of telescopes and satellites we 
can observe dark clouds at essentially all wavelengths possible, ranging from 7-rays 
to radio frequencies. Especially useful for studying star-forming regions is radio, 
submillimeter, and far-infrared (IR) emission. We have learned that all star forma- 
tion occurring in the Milky Way is associated with these dark clouds of molecular 
hydrogen and dust. These regions are sufficiently dense and well-shielded against 
the dissociating effects of interstellar ultraviolet radiation so that hydrogen atoms 
bind together to form molecules. Molecular hydrogen is a homonuclear molecule, 
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as a consequence its dipole moment vanishes and it radiates extremely weakly. 
Direct detection of cold interstellar H2 therefore is generally possibly only through 
ultraviolet (UV) absorption studies. Due to atmospheric opacity these studies can 
only be done from space, and are limited to pencil-beam measurements of the 
absorption of light from bright stars or active galactic nuclei. Note that rotational 
and ro-vibrational emission lines from H2 have also been detected in the infrared, 
both in the Milky Way and in other galaxies. However this emission comes from 
gas that has been strongly heated by shocks or radiation, and it traces only a small 
fraction of the total H 2 mass (e.g. van der Werf 20001. Due to these limitations, 
the most common tool for studying the molecular ISM is radio and sub-millimeter 
emission either from dust grains or from other molecules that tend to be found in 
the same locations as H2. The most important tracer molecule is CO, but also 
HCN, N2H+, NH3 and others are used depending on the density and temperature 
range of interest. 

Location and mass growth of young stars are intimately coupled to the dy- 
namical properties of their parental clouds. Stars form by gravitational collapse of 
shock-compressed density fluctuations generated from the supersonic turbulence 



ubiquitously observed in molecular clouds (e.g. Mac Low and Klessen 2004 Mc 



Kee and Ostriker 20071. Once a gas clump becomes gravitationally unstable, it 



begins to collapse and the central density increases considerably, giving birth to a 
protostar. In this dynamic picture, star formation takes place roughly on a free- 
fall timescale, as opposed to the "standard" model of the inside-out collapse of 
singular isothermal spheres, where core formation is dominated by the ambipolar 



diffusion timescale (Shu et al. 19871. Supersonic turbulence creates a highly tran- 



sient and inhomegeneous molecular cloud structure which is characterized by large 
density contrasts. Some the high-density fluctuations exceed the critical mass for 
gravitational contraction. The collapse of these Jeans-unstable cores leads to the 
formation of individual stars and star clusters. In this phase, a nascent proto- 
star grows in mass via accretion from the infalling envelope until the available gas 
reservoir is exhausted or stellar feedback effects become important and remove the 
parental cocoon — a new star is born. 

The structure of these lecture notes is as follows. First, we summarize the 
global properties of molecular clouds in Section [2] then we introduce basic aspects 
of interstellar turbulence in Section [3j and turn to the small-scale characteristics 
of molecular clouds in Section |4j where we discuss the properties of the cloud 
cores that are the direct progenitors of individual stars and binary systems. The 
statistical properties of stars and star clusters are the focus of Section 5l and we 
end with a critical analysis of different star formation theories in Section 6| By and 
large, these notes are based on two review articles, Mac Low and Klessen (2004) 



and Klessen et al. (2011) 



2 Global Properties of Molecular Clouds 



Molecular clouds in the Milky Way display a number of common properties. First, 
when studied with high spatial resolution clouds, they exhibit extremely complex 
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Fig. 1. The Cygnus X region is the richest and most massive complex of high-mass star 



formation at a distance lower than 3kpc. Schneider et al. (2006 20071 showed that the 



Cygnus X region constitutes a large scale network of GMCs at a common distance of 
~ 1.7 kpc. Signposts of recent and ongoing (high-mass) star formation (H II regions, 
outflow activity, IR-sources such as S106 IR, DR21, W75N, and GL2591) are ubiquitous. 



See Schneider et al. ( 2006 1 for a review on Cygnus X. The region also contains several 



OB clusters (Cyg OB14, 6, and 8) including the richest known OB cluster of the Galaxy, 



Cyg OB2 (Knodlseder 20001. Image from Schneider et al. '2011) 



and often filamentary structure, with column densities and corresponding 3-D 
densities that vary by many orders of magnitude. See Figure [l] for an image of 
Cygnus X, one of the richest and most massive complex of high-mass star formation 
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in the Galaxy. See also Table [T] for appropriate numbers. We discuss the small- 
scale structure and its relation to star formation in Section [4] below. On large 
scales, howerver, i.e. when observed at great distance or with low resolution, to 
within factors of a few all molecular clouds seem to have a similar mean surface 
density of ~ 100 M Q pc -2 corresponding to 0.035 g cm -2 (Heyer et al. 2008 



Bolatto et al. 2008). The constant surface density of molecular clouds is known 



as one of the Larson ( 1981 ) relations. We note, however, that there are a number 



of caveats with these relations and their interpretation ( Ballcsteros-Paredes and 



Mac Low 2002). Second, the clouds all display linewidths much greater than 



expected from thermal motion at typical temperatures of 10 — 20 K. On average, 
the observed linewidth is related to the size of the cloud by 



<7iD = 0.5 



1.0 pc 



0.5 



km s 



(2.1) 



where ctid is the one-dimensional cloud velocity dispersion and L is the cloud 
size (|Solomon et al.[ |1987| |Heyer and Brunt| |2004| |Bolatto et aL| |2008[). This is 



another one of Larson's relations. The large non-thermal linewidths have been 
interpreted as indicating the presence of supersonic turbulence, since both the low 
observed star formation rate (see below) and the absence of inverse P-Cygni line 
profiles indicates that they are not due to large-scale collapse. If one adopts this 
interpretation, then from these two observed relations one can directly deduce 
the third of Larson's relations, which is that giant molecular clouds have virial 



parameters (Bcrtoldi and McKee 1992) 



vir " GM ' 



(2.2) 



where M is the cloud mass and G is the gravitational constant. This indicates 
that these clouds are marginally gravitationally bound, but with enough internal 
turbulence to at least temporarily prevent global collapse. 

The presence of supersonic turbulence in approximate virial balance with self- 
gravity indicates that in molecular clouds the turbulent and gravitational energy 
densities are of the same order of magnitude, and both generally exceed the thermal 
energy density by large factors. If molecular clouds form in large-scale convergent 
flows, as we argue below, then surface terms from ram pressure can also be sig- 



nificant and need to be considered in the virial equations ( Ballesteros-Paredes 



2006). Also magnetic fields are important for the overall energy balance. The gas 



in molecular clouds is a weakly ionized plasma that is tied to magnetic field lines. 



Observations using Zeeman splitting (Crutcher 1999 Troland and Crutcher 2008 



Crutcher et al. 2010a) and the Chandrasekhar-Fermi effect (Lai et al. 2001 2002 



indicate that the field strength lies in the range from a few to a few tens of /iG. 
The exact values vary from region to region, but in general the magnetic energy 
density appears comparable to the gravitational and turbulent energy densities as 
expected from energy equipartition arguments. The dynamical importance of the 
field can be expressed in terms of the magnetic criticality. If the magnetic field 
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Table 1. Physical properties of molecular cloud and cores" 



molecular cluster- protostellar 
cloud forming cores 

clumps 



Size (pc) 


2-20 


0.1 - 


- 2 


<0.1 


Mean density (H2 cm" 3 ) 


10 2 - 10 3 


10 3 


- 10 5 


> 10 5 


Mass (M Q ) 


10 2 - 10 6 


10- 


- 10 3 


0.1 - 10 


Temperature (K) 


10-30 


10- 


- 20 


7-12 


Line width (kms" 1 ) 


1 - 10 


0.5- 


-3 


0.2-0.5 


RMS Mach number 


5-50 


2- 


15 


0-2 


Column density 










(g cm" 2 ) 


0.03 


0.03 


- 1.0 


0.3 - 3 


Crossing time (Myr) 


2-10 


<1 




0.1 - 0.5 


Free-fall time (Myr) 


0.3-3 


0.1 - 


- 1 


<0.1 


Examples 


Taurus, 


L1641, 


B68, L1544 




Ophiuchus 


L1709 





a Adapted from |Cernicharo| | |1991| l and |Bergin and Tafalla| | |2007[ r 



threading a cloud is sufficiently strong, then it cannot undergo gravitational col- 
lapse no matter what external pressure is applied to it, as long as it is governed by 
ideal magnetohydrodynamics (MHD). A cloud in this state is called subcritical. In 
contrast, a weaker magnetic field can delay collapse, but can never prevent it, and 
a cloud with such a weak field is called supercritical (IMouschovias and Spitzerl 



1976). Observations indicate the molecular clouds (Padoan and Nordlund 1999) 



and cores (Crutcher et al. 2009, 201 0a|b ), i.e. the high-density regions within the 
clouds that turn into individual stars or small stellar aggregates, are usually su- 
percritical with values within a factor of a few of the critical one (see also Section 
421. 



3 Interstellar Turbulence 

At this point, we need to digress from our discussion of molecular cloud prop- 
erties and turn our attention to the statistical description of turbulent flows. In 
particular, we focus on the differences between supersonic, compressible (and mag- 
netized) turbulence, which is characteristic of the interstellar medium (ISM), and 
the more commonly studied incompressible turbulence, which describes terrestrial 
flows, such as the motion of air in the Earth's atmosphere or the flow of water in 
rivers and oceans. For the purpose of our discussion, it is sufficient to think of tur- 
bulence as the gas flow resulting from random motions at many scales, consistent 
with the simple scaling relations discussed above. For a more detailed discussion of 
the complex statistical characteristics of turbulence in general, we refer the reader 
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to the book by Lesieur (1997), or for a thorough account of ISM turbulence to the 



reviews by Elmegreen and Scalo (2004) and Scalo and Elmegreen (2004). 



3.1 Subsonic, Incompressible Turbulence 

So far, most theoretical studies of turbulence treat incompressible turbulence. 
Root-mean-square (rms) velocities are subsonic, and the density remains almost 
constant. Dissipation of energy occurs primarily in the smallest vortices, where 



the dynamical scale I is shorter than the length on which viscosity acts v^. Kol 



mogorov ( 1941 ) described a heuristic theory based on dimensional analysis that 
captures the basic behavior of incompressible turbulence surprisingly well, al- 
though subsequent work has refined the details substantially. He assumed tur- 
bulence driven on a large scale L, forming eddies at that scale. These eddies 
interact to form slightly smaller eddies, transferring some of their energy to the 
smaller scale. The smaller eddies in turn form even smaller ones, until energy has 
cascaded all the way down to the dissipation scale ?7k- 

In order to maintain a steady state, equal amounts of energy must be trans- 
ferred from each scale in the cascade to the next, and eventually get dissipated, 
at a rate 

E = r]v 3 /L, (3.1) 

where 77 is a constant determined empirically. This leads to a power-law distri- 
bution of kinetic energy E oc v 2 oc fc~ n / 3 , where k = 2n/£ is the wavenumber, 
and density does not enter because of the assumption of incompressibility. The 
resulting differential energy spectrum is E(k)dk oc k~ 5 / 3 dk. Most of the energy 
remains near the driving scale, while energy drops off steeply below tjk- Because 
of the apparently local nature of the cascade in wavenumber space, the viscosity 
only determines the behavior of the energy distribution at the bottom of the cas- 
cade below ?7k, while the driving only determines the behavior near the top of the 
cascade at and above L. The region in between is known as the inertial range, 
in which energy gets transported from one scale to the next without influence 
from driving or viscosity. The behavior of the flow in the inertial range can be 
studied regardless of the actual scale at which L and tjk he, so long as they are 
well separated. Certain statistical properties of incompressible turbulence, such 
as structure functions S p (r) = ({v(x) — v(x + r)} p ) for example, have been suc- 
cessfully modeled by assuming that dissipation occurs in filamentary vortex tubes 



( She and Leveque 1994[ ) . 

Gas flows in the ISM, however, vary from this idealized picture in three impor- 
tant ways. First, they are highly compressible, with Mach numbers M. ranging 
from order unity in the warm, diffuse ISM, up to as high as 50 in cold and dense 
molecular clouds. Second, the equation of state of the gas is very soft due to 
radiative cooling, so that pressure P oc p 1 with the polytropic index falling in the 



range 0.4 < 7 < 1.2 as a function of density and temperature (e.g., Scalo et al. 



1998 Ballesteros-Paredes et al. 1999 Spaans and Silk 2000). Third, the driving 



of the turbulence is not uniform, but rather it is inhomogeneous and comes from 
a variety of sources that act on a range of different scales. 
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3.2 Supersonic, Compressible Turbulence 



Supersonic flows in highly compressible gas create strong density perturbations. 



Early attempts to understand turbulence in the ISM (e.g. von Weizsacker 1943 



von Weizsacker 1951 ) were based on insights drawn from incompressible turbu- 



lence. An attempt to analytically derive the density spectrum and resulting grav- 
itational collapse criterion was first made by |Chandrasekhar| ( |1951ajb| . This work 
was followed up by several authors, culminating in work by Sasao (1973) on den- 
sity fluctuations in self-gravitating media whose interest has only been appreciated 



recently. Larson (1981) qualitatively applied the basic idea of density fluctuations 



driven by supersonic turbulence to the problem of star formation. Bonazzola et al.] 
( 1973 ) used a renormalization group technique to examine how the slope of the 



turbulent velocity spectrum could influence gravitational collapse. This approach 
was combined with low-resolution numerical models to derive an effective adiabatic 
index for subsonic compressible turbulence by Panis and Perault (1998). Adding 
to the complexity of the problem, the strong density inhomogeneities observed in 
the ISM can be caused not only by compressible turbulence, but also by thermal 



phase transitions ( Field et al.| 1969 McKee and Ostriker[ 1977, Wolfire et al. 



1995), or gravitational collapse (see, e.g., Kim and Ostriker| 2001 for the Galaxy, 



Schleicher et~aT| |2010l |Sur et al.[ |2010[ |Federrath et al.[ \20U\ for applications 



or 

to high redshift halos). 

In supersonic turbulence, shock waves offer additional possibilities for dissi- 
pation. Shock waves can also transfer energy between widely separated scales, 
removing the local nature of the turbulent cascade typical of incompressible tur- 
bulence. The spectrum may shift only slightly, however, as the Fourier transform 
of a step function representative of a perfect shock wave is k~ 2 . For purely shock- 
dominated turbulence, the so called Burgers' turbulence, the resulting spectrum 
is E(k)dk oc k~ 2 dk. Here, certain statistical characteristics can be understood 
by assuming dissipating occurs in sheetlike shocks ( |Boldyre vi 



2002 Padoan et al. 



2004 Schmidt et al. 2008). Numerical simulations, as well as observations in the 



interstellar medium, indicate an spectral slope that lies somewhere in between 
Kolmogorov's value of —5/3 and the value of —2 for Burgers' turbulence (see, e.g., 
Federrath et all 120091 12010 1. 



3.3 Origin of ISM Turbulence 



The physical origin of turbulence in the ISM is not fully understood yet. In 
particular, the question as to whether it is injected from outside or driven by 
internal sources, i.e. whether it is driven on large or small scales, is still subject 
to considerable debate. We favor the first assumption as observations indicate 
that molecular cloud turbulence is always dominated by the largest-scale modes 



accessible to the telescope ( Ossenkopf and Mac Low 2002 Brunt et al. 2009 ) . 
In addition, the amount of turbulence in molecular clouds with no or extremely 
low star formation like the Magdalena cloud or the Pipe nebula, is significant 
and broadly comparable to the level of turbulence observed towards star forming 
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clouds. Both facts seem difficult to reconcile with turbulence being driven from 
internal stellar sources. Instead, we argue that it is the very process of cloud 
formation that drives its internal motions by setting up a turbulent cascade that 
transports kinetic energy from large to small scales in a universal and self-similar 
fashion. 

The driving sources for turbulence on the largest scales can be diverse, rang- 
ing from the accretion of gas of extragalactic origin (Klessen and Hennebelle[ 



2010) to the occurance of convergent flows of atomic gas triggered by spiral den- 



sity waves (Walder and Folini 
Heitsch et all |2006bl 120081 iHennebelle et al 



1996 1998 Vazquez-Semadeni et al. 2007 2009 



nova ex plosions (|Mac Low and Klessen 2004 



regions (Matzner 2002 Krumholz et al. 2006 Peters et al. 



2008 Banerjee et al. 2009), super- 



Dib et al. 2006), or expanding HII 



2008 Gritschneder 



et al. 2009 ). Here it is the very process of cloud formation that drives the internal 
turbulence. Some models have also investigated molecular cloud turbulence that 



is driven on small scales by internal sources such as stellar winds and outflows ( Li 



and Nakamura , 2006 Nakamura and Li 2007 Wang et al. 2010 ) but it is unlikely 



that these have a significant effect on the largest scales within clouds (Mac Low 



and Klessen 2004 Banerjee et al. 2007 Brunt et al. 20091. It is an appealing 



hypothesis that molecular clouds form at the stagnation points of large-scale con- 
vergent flows ( Ballesteros-Paredes et al. 1999 2005). As the density goes up the 



gas can cool efficiently, turn from being mostly atomic to molecular, and shield 
itself from the external radiation field. As long as the convergent flow continues 
to deliver fresh material the cloud grows in mass and is confined by ram pressure. 
Recent numerical simulations attempting to form molecular clouds from diffuse 



gas ( 


Vazquez-Semadeni et al. 2007 


Banerjee et al. 


2009 


) by considerin; 



Heitsch et al. 2008 Hennebellc et al. 2008 



have indeed shown that this process is sufficient to sustain a substantial degree of 
turbulence in the forming cloud. 



3.4 Relation to Star Formation 

Altogether, we propose the following picture. On large scales, the turbulence in 
molecular clouds is highly supersonic. We know that the density contrast for 
isothermal gas scales with the Mach number M. to the second power, that means 
with M. « 10 we expect density contrasts of roughy 100. This is clearly observed 
in molecular clouds, where the mean density is around 100 particles per cubic cen- 
timeter and where the high-density cores exceed values of 10 4 cm -3 and more (see 
Table [l]). When zooming in on cluster-forming cloud cores (or their not-yet-star- 
forming counterparts, the so-called infrared dark clouds) one still observes RMS 
Mach numbers of ~ 5, which still leads to localized density fluctuations of a factor 
of 25 on average. As indicated earlier, some of these fluctuations may exceed the 
critical mass for gravitational collapse to set in. The presence of turbulence, there- 
fore, makes the cluster-forming core to break apart into smaller units. It fragments 
to build up a cluster of stars with a wide range of masses rather than forming one 
single high-mass star. We call this process gravoturbulent fragmentation, because 
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Fig. 2. Simple cartoon picture of the turbulent energy spectrum, i.e. of the kinetic 
energy carried by modes of different wave numbers k, and their relation to different 
cloud structures (see also Table [TJ. Turbulence is driven on large scales comparable to 
the size L of the cloud and is dissipated on very small scales r/K- 



turbulence generates the distribution of clumps initially and then gravity selects 
a subset of them for star formation. For a more detailed account, see |Mac Low| 
and Klessen (2004). Finally, when focussing on low-mass cores the velocity field 

(e.g. |Goodman et al.[ |1998[ |Bergin and Tafalla| |2007[) 



becomes more coherent 
and the turbulence subsonic 



This defines the sonic scale at around ~ 0.1 pc. 



Such structures are no longer subject to gravoturbulent fragmentation and are 
the direct progenitors of individual stars or binary systems. We note, however, 
that gravitational fragmentation may still occur within the protostellar accretion 
disk that builds up in the center of the core due to angular momentum conser- 
vation. This process is thought to produce close binaries (see, e.g., |Bodenheimer 



1995 Machida 2008). The fact that the observed velocity dispersion approaches 



the thermal value as one zooms in on smaller and smaller scales, as expressed in 
Larson's second relation (equation 2.1), is a direct consequence of the negative 
power-law slope of the turbulent energy spectrum, as illustrated in Figure [2] 



4 Molecular Cloud Cores 



In this section focus on the small-scale characteristics of molecular clouds and 
discuss the properties of the low-mass cores that are the immediate progenitors of 
individual stars. 
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4.1 Statistical Properties 



Emission line observations and dust extinction maps of molecular clouds reveal 
extremely complex morphological structure with clumps and filaments on all scales 
accessible by present day telescopes. Typical parameters of different regions in 
molecular clouds are listed in Table [T] The volume filling factor of dense clumps 
with density contrast n/(n), where n is the local density and (n) the average 
density of the cloud, is rather small. It is of order of a few per cent for densities 
10 5 cm" 3 dBlitzl 1 19931 iMcKeel Il999l IWilliams et all |2000l iBeuther et al. 



n > 



2000). As discussed above, this is a direct consequence of the fact that ISM 
turbulence is highly supersonic on large scales. It is important to note that star 
formation always occurs in the densest regions within a cloud, so only a small 
fraction of molecular cloud matter is actually involved in building up stars, while 
the bulk of the material remains at lower densities. This is the key to explaining 
the low star formation efficiencies in Galactic molecular clouds. 

The mass spectrum of clumps in molecular clouds appears to be well described 
by a power law, 

dN , s 

^ m < ^ 
with the exponent being in the range —1.3 < a < —1.8, indicating self-similarity 



( Stutzki and Gusten[|1990[|Williams et al.[|1994[|Kramer et al.[|1998[ ). There is no 
natural mass or size scale between the lower and upper limits of the observations. 
The smallest observed structures are protostellar cores with masses of a few solar 
masses or less and sizes of ^ 0.1 pc. Given the uncertainties in determining the 
slope, it appears reasonable to conclude that there is a universal mass spectrum 
for the clumps within a molecular cloud, and that the distribution is a power 
law within a mass range of three orders of magnitude, i.e. from 1 Mq to about 
1000 M . Hence, it appears plausible that the physical processes that determine 
the distribution of clump masses are rather similar from cloud to cloud, and that 
they are closely related to the universal nature of turbulent flows and thermal 
instability acting on self-gravitating gas. 

Most of the objects that enter in the above morphological analyses are not 



gravitationally bound ( 


Stutzki and Glisten 


et al. 


2005 


Dib et al. 


20071. It is interesting 



2005 Klessen 



as one probes smaller and smaller scales and more and more bound objects. When 
considering prestellar cores, which are thought to be the direct progenitors of 
individual stars or small multiple systems, then the mass function is well described 
by a double power law fit dN/dm cx mT a following a — 2.5 above ~ 0.5 M and 
a = 1.5 below. The first large study of this kind was published by |Motte et al.| 



( 1998 ), for a population of submillimetre cores in p Oph. Using data obtained with 
IRAM, they discovered a total of 58 starless clumps, ranging in mass from 0.05 M e 



to ~ 3 M . Similar results are obtained from the Serpens cloud ( |Testi and Sargent 



et al. 



1998), for Orion B North (Johnstone et al. 2001) and Orion B South (Johnstone 



2006 ), or for the Pipe Nebula (Lada et al. 2006 1. Currently all observational 



data (Motte et al. 1998 Testi and Sargent 1998 Johnstone et al. 2000 2001 
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2006 



2007 



Nutter and Ward-Thompson) |2007| |Alves et al.[ |2007[ |Di Francesco et al.[ 



Ward-Thompson et al. 



2007 



Lada et al. 2008 ) reveal that the mass function 



of prestellar cores is strikingly similar in shape to the stellar initial mass function, 
the IMF. To reach complete overlap one is required to introduce a mass scaling or 
efficiency factor in the range 2 to 10, which differs in different regions. An exciting 
interpretation of these observations is that we are witnessing the direct formation 
of the IMF via fragmentation of the parent cloud. However, we note that the 
observational data also indicate that a considerable fraction of the prestellar cores 
do not exceed the critical mass for gravitational collapse, much like the clumps on 
larger scale. The evidence for a one-to-one mapping between prestellar cores and 
the stellar mass thus is by no means conclusive (see, e.g. 



Clark et al. 2007) 



4.2 Individual Cores 



Density Structure. The density structure of prestellar cores is typically estimated 
through the analysis of dust emission or absorption using near-IR extinction map- 
ping of background starlight, millimeter/submillimeter dust continuum emission, 
or dust absorption against the bright mid-IR background emission (jBergin and| 



Tafalla 2007 ) . A main characteristic of the density profiles derived with the above 



techniques is that they require a central flattening. The density gradient of a core 
is flatter than r _1 within radii smaller than 2500 — 5000 AU, with typical central 
densities of 10 5 - 10 6 cixT 3 (iMotte et all Il998l IWard-Thompson et all Il999h. 



A popular approach is to describe these cores as truncated isothermal (Bonnor- 
Ebert) sphere (Ebert 1955 Bonnor 1956), that often (but not always) provides 



a good fit to the data (Bacmann et al. 2001 Alves et al| 2001 Kandori et al. 



2005). These are equilibrium solutions of self-gravitating gas spheres bounded by 
external pressure. However, such density structure is not unique. Numerical calcu- 
lations of the dynamical evolution of supersonically turbulent clouds show that the 
transient cores forming at the stagnation points of convergent flows exhibit similar 



morphology despite not being in dynamical equilibrium ( jBallesteros-Paredes et al 
20031). 



Thermal Stucture. The kinetic temperature of the dust and gas components in 
a core is regulated by the interplay between various heating and cooling processes. 
At high densities (> 10 5 cm~ 3 ) in the inner part of the cores, the gas and dust 



have to be coupled thermally via collisions (Goldsmith and Langer 1978 Burke 



and Hollenbach 1983 Goldsmith 20011. At lower densities, which correspond to 



the outer parts of the cores, the two temperatures are not necessary expected to 
be the same. Thus, the dust and gas temperature distributions need to be inferred 
from observations independently, large-scale studies of dust temperature show 
that the grains in starless cores are colder than in the surrounding lower-density 
medium. Far-IR observations toward the vicinity of a number of dense cores pro- 
vide evidence for flat or decreasing temperature gradients with cloud temperatures 
of 15 - 20 K and core values of 8 - 12 K (IWard-Thompson et all 120021 |T6th et al. 



2004). These observations are consistent with dust radiative transfer modeling 



in cores illuminated by interstellar radiation field ( jlanger et al. 2005 Keto and 
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Fig. 3. Formation and growth of molecular cloud cores by thermal instability triggered by 
a large-scale convergent flow. A small cold condensate grows from the thermally unstable 
warm neutral medium by outward propagation of its boundary layer. Coalescence and 
merging with nearby clumps further increases its mass and size. The global gravitational 
potential of the proto-cloud enhances the merging probability with time. The images 
show 2D slices of the density (logarithmic colour scale) and the gas velocity (indicated as 
arrows) in the plane perpendicular to the large scale flow. From Banerjee et al. \2009 ). 



Field, 2005 Stamatellos et al. 2007 1, where the dust temperature is ~ 7K in the 
core center and increases up to 16 K in the envelope. The gas temperature in 
molecular clouds and cores is commonly infered from the level excitation of simple 
molecules like CO and NH3 (Evans 1999 Walmsley and Ungerechts 19831. One 
finds gas temperatures of 10 — 15 K, with a possible increase toward the lower 
density gas near the cloud edges. It is believed that the gas heating in prestellar 
cores mostly occurs through ionization by cosmic rays, while the cooling is mainly 
due to line radiation from molecules, especially CO (Goldsmith and Langer 1978 1. 
Altogether, the fact that prestellar cores are cold and roughly isothermal with at 
most a modest increase in temperature from the center to the edge is consistent 



with numerical models of cores forming from thermal instability (Heitsch et al. 
2006a| |Keto and Caselli[ |2008| |Banerjee et"aTj |2009[ ), see also Figure |3J 



Chemical Stucture. Maps of integral line intensity can look very different for dif- 
ferent molecular tracers. In particular, the N2H 4 " and NH3 emission more closely 
follows the dust emission while the C 18 and CS emission appears as a "ring- 
like" structure around the dust emission maximum ( Bergin et al. 2002 Tafalla 



et al.[ |2002| |Lada et al.| |2003[ |Maret et alj |2007[ ). For illustration see Figure |4j 



The common theoretical interpretion of these data is that carbon-bearing species, 
represented by CO and CS freeze-out on the dust grains from the gas while the 
abundances of nitrogen- hydrogen bearing molecules, N 2 H + and NH 3 , either stay 
constant or decay more slowly. At the same time, chemical models of prestellar 
cores predict that molecules in the core envelope have to be destroyed by inter- 
stellar UV field (iPavlyuchenkov et all 120061 lAikawa et all 120081). The chemical 



stratification significantly complicates the interpretation of molecular line obser- 
vations and again requires the use of sophisticated chemical models which have to 
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Fig. 4. Maps of molecular line emission from C ls O, N2H + , and CS superimposed on 



a dust extinction maps of the dark cloud core Barnard 68 (Alves et al. 2001 Bergin 



|et al.| |2002| |Lada et al. 1 12003| ) . The three images illustrate the effects of depletion onto 
grains in the high-density central region of the core. N2H + is the least and CS the most 
depleted species. Image courtesy of E. A. Bergin. 



be coupled to the dynamical evolution (see, e.g., Glover and Mac Low 2007a|b 



Glover et al. 20101. From observational side, the freeze-out of many molecules 



makes it difficult to use their emission lines for probing the physical conditions in 
the inner regions of the cores. At the same time, the modeling of the chemical 
evolution can provide us with the important parameters of the cores. For example, 
the level of CS depletion can be used to constrain the age of the prestellar cores 
while the deficit of CS in the envelope can indicate the strength of the external UV 
field (Bergin and Tafalla 2007). In any case, any physical interpretation of the 



molecular lines in prestellar cores has to be based on chemical models and should 
do justice to the underlying density and velocity pattern of the gas. 

Kinematic Stucture. In contrast to the supersonic velocity fields observed in 
molecular clouds, dense cores have low internal velocities. Starless cores in clouds 
like Taurus, Perseus, and Ophiuchus systematically present spectra with close-to- 



thermal linewidths, even when observed at low angular resolution (Myers 1983 



Jijina et al. 1999). This indicates that the gas motions inside the cores are sub- 



sonic or at best transsonic, i.e. with Mach numbers < 2 (Kirk et al. 



2007 



Andre 



et al. , 2007 Rosolowsky et al. 2008 1 . In some cores also inward motions have been 



detected. They are inferred from the observation of optically thick, self-absorbed 
lines of species like CS, H2CO, or HCO + , in which low-excitation foreground gas 
absorbs part of the background emission. Typical inflow velocities are of order 
of 0.05 — O.lkm/s and are observed on scales of 0.05 — 0.15 pc, comparable to 
the observed size of the cores (Lee et al. 19991. The overall velocity structure 



of starless cores appears broadly consistent with the structure predicted by mod- 
els in which protostellar cores form at the stagnation points of convergent flows, 
but the agreement is not perfect. Simulations of core formation do correctly find 
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that most cores are at most transsonic (Klessen et al. 2005 Offner et al. 2008), 



but the distribution of velocity dispersions has a small tail of highly supersonic 
cores that is not observed. Clearly more theoretical and numerical work is needed. 
In particular, the comparison should be based on synthetic line emission maps, 
which requires to couple a chemical network and radiative transfer to the sim- 
ulated density profiles as discussed above. In addition, it is also plausible that 
the discrepancy occurs because the simulations do not include all the necessary 
physics such as radiative feedback and magnetic fields. Subsonic turbulence con- 
tributes less to the energy budget of the cloud than thermal pressure and so cannot 



provide sufficient support against gravitational collapse (Myers 1983 Goodman 



et al. 1998 Tafalla et al. 2006). If cores are longer lasting entities there must 



be other mechanisms to provide stability. Obvious candidates are magnetic fields 



( Shu et al. 


1987 


). However, they are usually not strong enough to provide suffi- 


cient support ( 


Crutcher et al. 


1999 


Crutcher and Troland 


2000 


Bourke et al. 


2001 


Crutcher et al. 


2009 


2010b 1 as discussed below. Most observed cores are 



thus likely to be evolving transient objects that never reach any equilibrium state. 
Magnetic Field Structure. Magnetic fields are ubiquitously observed in the 



interstellar gas on all scales (Crutcher et al. 2003 Heiles and Troland 2005). 



However, their importance for star formation and for the morphology and evolu- 
tion of molecular cloud cores remains controversial. A crucial parameter in this 
debate is the ratio between core mass and magnetic flux. In supercritical cores, 
this ratio exceeds a critical value and collapse can proceed. In subcritical ones, 
magnetic fields provide stability (Spitzer 1978 Mouschovias 1991a|b ). Measure- 
ments of the Zeeman splitting of molecular lines in nearby cloud cores indicate 
mass-to-flux ratios that lie above the critical value, in some cases only by a small 



margin but very often by factors of many if non-detections are included ( Crutcher 



1999 Bourke et al. 2001 Crutcher et al. 2009 2010a). The polarization of dust 



emission offers an alternative pathway to studying the magnetic field structure 
of molecular cloud cores. MHD simulations of turbulent clouds predict degrees of 
polarization between 1 and 10%, regardless of whether turbulent energy dominates 



over the magnetic energy (i.e. the turbulence is super- Alfvenic) or not (Padoan 



and Nordlund 1999 


Padoan et al. 




2001 


. However, converting 


magnetic field strength is very difficult 


Heitsch et al. 




2001b) 



Altogether, the 

current observational finding imply that magnetic fields must be considered when 
studying stellar birth, but also that they are not the dominant agent that deter- 
mines when and where stars form within a cloud. Magnetic fields appear too weak 
to prevent gravitational collapse to occur. 

This conclusion means that in many cases and to reasonable approximation 
purely hydrodynamic calculations are sufficient for star formation simulations. 
However, when more precise and quantitative predictions are desired, e.g. when 
attempting to predict star formation timescales or binary properties, it is necessary 
to perform magnetohydrodynamic (MHD) simulations or even consider non-ideal 
MHD. The latter means to take ambipolar diffusion (drift between charged and 
neutral particles) or Ohmic dissipation into account. Recent numerical simula- 
tions have shown that even a weak magnetic field can have noticeable dynamical 



Ralf S. Klessen: Star Formation in Molecular Clouds 



15 



50 



|-50 



-100 



100 



50 



R.A. Offset {") 



-50 



-100 



Fig. 5. Polarization vector map of the central region of the Serpens cloud core, super- 
posed on the total intensity images in logarithmic scaling. The area of the image is 220" 
x 220". From\Sugitani et al.\ (Mlffy. 



effects. It can alter how cores fragment (Price and Bate 2007b 2008| Hennebelle 


and Fromang 2008 


Hennebelle and Teyssier[ |2008 , 


Hennebelle and Ciardi| 2009 , 


Hennebelle et al. 1 


2011 Peters et al. 


201ip , change the coupling between stel- 



lar feedback processes and their parent clouds (Nakamura and Li 2007 Krumholz 



et al. 2007b ), influence the properties of protostellar disks due to magnetic braking 
(Price and Bate 2007a Mellon and Li 2008), or slow down the overall evolution 



(Heitsch et al. 2001a 



5 Statistical Properties of Stars and Star Clusters 



In order to better understand how gas turns into stars, we also need to discuss 
here some of the key properties of young stellar systems. We restrict ourself to a 
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discussion of the star formation timescale, the spatial distribution of young stars, 
and the stellar initial mass function (IMF). We note, however, that other statistical 
characteristics, such as the binary fraction, their relation to the stellar mass, and 
the orbital parameters of binary stars are equally important for distinguishing 
between different star formation models. 



5.1 Star Formation Timescales 



The star formation process in molecular clouds appears to be fast (Hartmann et al. 



2001 Elmegreen 20071. Once the collapse of a cloud region sets in, it rapidly forms 
an entire cluster of stars within 10 6 years or less. This is indicated by the young 
stars associated with star forming regions, typically T Tauri stars with ages less 
10 6 years (Gomez et al. 1992 Greene and Meyer 1995 Carpenter et al. 



than 



1997), and by the small age spread in more evolved stellar clusters (Hillenbrand 



1997 



Palla and Stahler 1999). Star clusters in the Milky Way also exhibit an 



amazing degree of chemical homogeneity (in the case of the Pleiades, see Wilden 



et al. , 2002 ), implying that the gas out of which these stars formed must have been 



chemically well-mixed initially (see also de Avillez and Mac Low 
and Linl 120031). 



2002 Klessen 



5.2 Spatial Distribution 

The advent of sensitive infrared detectors in the last decade has made it possible 
to perform wide-area surveys. These have led us to recognize that most stars form 
in clusters and aggregates of various size and mass scales, and that isolated or 



widely distributed star formation is the exception rather than the rule (Lada and 
Lada 2003 1 . The complex hierarchical structure of molecular clouds (Figure |lj) 
provides a natural explanation for this finding. 

Star-forming molecular cloud cores vary enormously in size and mass. In small, 
low-density, clouds stars form with low efficiency, more or less in isolation or scat- 
tered around in small groups of up to a few dozen members. Denser and more 
massive clouds may build up stars in associations and clusters of a few hundred 
members. This appears to be the most common mode of star formation in the 
solar neighborhood (Adams and Myers 2001). Examples of star formation in 



small groups and associations are found in the Taurus- Aurigae molecular cloud 



( Hartmann 


2002 


). Young stellar { 


Chamaeleon I ( 


Persi et al. 2000) 



Each of these clouds is at a distance of about 130 to 160 pc from the Sun. Like 
most of the nearby young star forming regions they appear to be associated with 



a ring- like structure in the Galactic disk called Gould's belt (Poppel 1997). 



The formation of dense rich clusters with thousands of stars is rare. The closest 



region where this happens is the Orion Nebula Cluster in L1641 (Hillenbrand 1997 



Hillenbrand and Hartmann 1998), which lies at a distance of 410 pc (Sandstrom 



etal. 2007 Mentenetal. 2007 Hirotaetal. 2007 Caballero 2008). A rich cluster 
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Fig. 6. Comparison of clusters of different masses scaled to same relative distance. The 
cluster in the upper left corner is the Orion Nebula Cluster ( McCaughrean 2001 1 and the 



one at the lower left is NGC 3603 (Brandl et al. 1999 1, both observed with the Very Large 
Telescope at infrared wavelength. The large cluster in the center is 30 Doradus in the 
LMC observed with the Hubble Space Telescope (courtesy of M. J. McCaughrean). The 
total mass increases roughly by a factor of ten from one cluster to the other. Composite 
image courtesy of H. Zinnecker. 



somewhat further away is associated with the Monoceros R2 cloud (Carpenter 
et al. 1997) at a distance of ~ 830 pc. The cluster NGC 3603 is roughly ten 
times more massive than the Orion Nebula Cluster. It lies in the Carina region, 
at about 7kpc distance. It contains about a dozen O stars, and is the nearest 



object analogous to a starburst knot (Brandl et al. 1999 Moffat et al. 2002) 



To find star-forming regions building up hundreds of O stars one has to look 
towards giant extragalactic Hn regions, the nearest of which is 30 Doradus in 
the Large Magellanic Cloud, a satellite galaxy of our Milky Way at a distance at 
55kpc. The giant star forming region 30 Doradus is thought to contain up to a 
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Fig. 7. Stellar mass spec- 
trum in different nearby clusters 
(black symbols: Orion Nebula 
Cluster, green: Pleiades, blue: 
M35) and its description by a 
three-component power law (red 
lines with overall uncertainties 
indicated by the hatched re- 
gion). From \Kroupa\ 



log 10 m [MJ 



hundred thousand young stars, including more than 400 O stars (Hunter et al. 



1995 


Walborn and Blades 


1997 


Townsley et al. 2006 



magnitude in scale, ranging from isolated single stars to massive young clusters 
with several 10 4 stars. This enormous variety of star forming regions appears to 
be controlled by the competition between self-gravity and opposing agents such 
as the turbulence in the parental gas clouds, its gas pressure and magnetic field 
content. 



5.3 The Stellar Initial Mass Function 

Mass is the most important parameter determining the evolution of individual 
stars. Massive stars with high pressures at their centers have strong nuclear fusion, 
making them short-lived but very luminous, while low-mass stars are long-lived 
but extremely faint. For example, a star with 5 M only lives for 2.5 x 10 7 yr, while 
a star with 0.2 M survives for 1.2 x 10 13 yr, orders of magnitude longer than the 
current age of the universe. For comparison the Sun with an age of 4.5 x 10 9 yr 
has reached approximately half of its life span. The relationship between mass and 



luminosity is quite steep with roughly L cx M 3 2 (Kippenhahn and Weigert 1994). 
During its short life a 5 M star will shine with a luminosity of 1.5 x 10 4 L , while 
the luminosity of an 0.2 M star is only ~ 10 -3 L . For reference, the luminosity 
of the Sun is 1L = 3.85 x lO^ergs -1 . 
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Explaining the distribution of stellar masses at birth, the so-called initial mass 
function (IMF), is a key prerequisite to any theory of star formation. The IMF 
has three properties that appear to be relatively robust in diverse environments, 
see Figure |7| These are the power law behavior dN/dm ex m a with slope a 2.3 
for masses m above about 1 M , originally measured by Salpeter ( 1955 ) , the lower 



mass limit for the power law and the broad plateau below it before the brown 

and the maximum mass of 



dwarf regime (Miller and Scalo 1979 Scalo 



stars at around 100 M© (Weidner and Kroupa 



1986) 



2004 



2006 Oey and Clarke 



Comprehensive reviews of cluster and field IMFs may be found in Scalo (1986), 



Kroupa (2002), Chabrier (2003), and Bastian et al. (2010) 



2005) 



At the extreme ends of the stellar mass spectrum, however, our knowledge 
of the IMF is limited. Because massive stars are very rare and short lived, only 
very few are sufficiently near to study them in detail and with very high spatial 
resolution, for example to determine multiplicity ( Zinnecker and Yorke 2007 ) . 



Low-mass stars and brown dwarfs, on the other hand, are faint, so they too are 
difficult to study in detail (Burrows et al. 2001). Such studies, however, are 
in great demand, because secondary indicators such as the fraction of binaries 
and higher-order multiples as function of mass, or the distribution disks around 
very young stars and possible signatures of accretion during their formation are 
probably better suited to distinguish between different star-formation models than 
just looking at the IMF. In contrast to the observational agreement on the IMF, 
at least above the substellar regime, there is still considerable disagreement on the 
theoretical side. The origin of the IMF is a major topic of theoretical research and 
we can only provide some basic arguments here. For more extended discussions, 
seelMac Low and Klessenl (|2004|), iBonnell et al.l (|2007|), iLarsonl (|2007|), orlMcKee 



and Ostriker (2007). 



6 Star Formation 



We begin this section with a brief discussion about the possible origin of low-mass 
stars, which was traditionally the focus on theoretical models of stellar birth. We 
argue that in the past years, we have seen a paradigm shift away from a slow, 
quasi-static picture of star formation to a faster and more dynamic approach that 



is based on the interplay between gravity and turbulence (see also Mac Low and 
Klessen, 2004). This provides a natural framework for a more consistent under- 
standing of star and star cluster formation. Next we give a critical account of 
the various proposals in the literature to explain the physical origin of the IMF 
within the context of the current dynamical theory of star formation. Finally, 
we turn our attention to the high- mass end of the IMF. Massive stars influence 
the surrounding universe far out of proportion to their numbers through ionizing 
radiation, supernova explosions, and heavy element production. Their formation 
requires the collapse of massive interstellar gas clouds with accretion rates exceed- 
ing 10" 



Mq yr" 1 (Beuther et al. 2002 Beltran et al 



masses before exhausting their nuclear fuel (|Keto and Wood 2006|) 
sons 



2006 ) to reach their final 
For these rea- 



the formation of massive stars is less well understood than the formation of 
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low mass stars. However, due to recent advances in algorithmic development and 
computational capabilities, which allow us to perform 3-dimensional, radiation- 
magnctohydrodynamic calculations of the collapse and subsequent star formation 



of massive cloud cores (see, e.g., Klessen et al. 2011 1, our understanding of the ori 



gin of the extreme ends of the stellar mass distribution has increased enormously. 
These developments are discussed in the last part of this Section. 



6.1 Dynamic Star Formation Theory 



The past decade has seen a paradigm shift in low-mass star formation theory. The 
general believe since the 1980's was that cores in low-mass star-forming regions 



evolve quasi-statically in magnetically subcritical clouds 1 
picture, gravitational contraction is mediated by ambipola 


Shu et al. 
r diffusion 


1987]. In this 
( Mouschovias 


1976 


1979 


Mouschovias and Paleologou 1981 


Mouschovias 


1991a 


I causing a re- 



distribution of magnetic flux until the inner regions of the core become supercritical 
and go into dynamical collapse. This process was originally thought to be slow, 
because in highly subcritical clouds the ambipolar diffusion timescale is about 
10 times larger than the dynamical one. However for cores close to the critical 
value, as is suggested by observations, both timescales are comparable. Numerical 
simulations furthermore indicate that the ambipolar diffusion timescale becomes 
significantly shorter for turbulent velocities similar to the values observed in nearby 



mura 



(Fatuzzo and Adams 


2002 


Heitsch et al. 


2004 


Li and Naka- 



2004 1 . The fact that ambipolar diffusion may not be a slow process under 



realistic cloud conditions, as well as the fact that most cloud cores are magnet- 



ically supercritical (Crutcher et al. 1999 Crutcher and Troland 2000 Bourke 



et al. 2001 Crutcher et al. 2009 ) has cast significant doubts on any magnetically- 
dominated quasi-static models of stellar birth. For a more detailed account on the 
shortcomings of the quasi-static, magnetically dominated star formation model, 
see 



Mac Low and Klessen ( 2004 \ 



For this reason, star-formation research has turned into considering supersonic 
turbulence as being on of the primary physical agents regulating stellar birth. 
The presence of turbulence, in particular of supersonic turbulence, has important 
consequences for molecular cloud evolution. On large scales it can support clouds 
against contraction, while on small scales it can provoke localized collapse. Turbu- 
lence establishes a complex network of interacting shocks, where dense cores form 
at the stagnation points of convergent flows. The density can be large enough 
for gravitational collapse to set in. However, the fluctuations in turbulent veloc- 
ity fields are highly transient. The random flow that creates local density en- 
hancements can disperse them again. For local collapse to actually result in the 
formation of stars, high density fluctuations must collapse on timescales shorter 
than the typical time interval between two successive shock passages. Only then 
are they able to 'decouple' from the ambient flow and survive subsequent shock 
interactions. The shorter the time between shock passages, the less likely these 
fluctuations are to survive. Hence, the timescale and efficiency of protostellar core 
formation depend strongly on the wavelength and strength of the driving source 
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turbulence creates a hierarchy of clumps while the whole region contracts, individual individual clumps collapse to form stars 

clumps collapse to form stars 








- 







in dense dusters clumps may merge while in dense dusters competitive mass growth 
collapsing -* contain multiple protostars becomes important 







in dense dusters N-body effects influence 
mass growth 






low-mass objects may become ejected 
-+ accretion stops 




feedback terminates star formation 







■•a 



result: star cluster, possibly with HII region 



Fig. 8. Star cluster formation in a turbulent molecular cloud core. 



(Klessen et al. 2000 Heitsch et al 


2001a Vazquez- Semadeni et al.| 


2003 Mac 


Low and Klessen 2004 Krumholz and McKee[ 2005[ Ballesteros-Paredes et al.[ 


2007 


McKee and Ostriker 


2007 1 , and accretion histories of individual protostars 


are strongly time varying ( 


Klessen 


2001b 


Schmeja and Klessen 


2004 


1. 



Altogether, we propose an evolutionary sequence as outlined in Figure [8} Fol- 
lowing the discussion in Section |3.4| star cluster formation takes place in massive 
cloud cores of several 10 2 to 10 3 solar masses with sizes of a few parsec and a veloc- 
ity dispersion of ~ 1 kins -1 (see also Table [TJ. In order to form a bound cluster, 
the potential energy must dominate the energy budget, meaning that the entire 
region is contracting. This is a typical feature of the structures formed by large- 
scale convergent flows (see Section 3.3). On the scales of a cluster-forming core, 
the internal turbulent motions are still supersonic with Mach numbers A4 ~ 5, 
and as a consequence there is a high degree of internal substructure with large 
density contrasts. Some of these density fluctuations are gravitationally unsta- 
ble in their own right and begin to collapse on timescales much shorter than the 
global contraction time. Recall that the free-fall time Tg scales with the density p 



as Tg oc p 



-1/2 
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Typically, the most massive fluctuations have the highest density and form a 
protostar in their center first. This nascent star can accrete from the immediate 
environment, but because it is located in a minimum of the cloud core's gravita- 
tional potential more gas flows towards it, and it can maintain a high accretion 
rate for a long time. In contrast, stars that form in lower-mass gas clumps typically 
can only accrete material from their immediate surrounding and not much beyond 



that (see, e.g., Klessen and Burkert 2000 2001 Klessen 2001 Bonnell et al. 



2004). Because this preferentially happens in the cluster outskirts, these processes 



naturally lead to mass segregation as we often observe in young clusters (see, e.g. 
[Hillenbrand] [19971 |Hillenbrand and Hartmamr} [T9981 for the Orion Nebula Clus- 
ter). In very dense clusters, there is the possibility that clumps merge while still 
accreting onto their central protostars. These protostars now compete for further 
mass growth from a common gas reservoir. This is the essence of the competitive 
accretion model (Bonnell and Bate 2006). In its original flavor, it was thought 



that this leads to the run-away growth of the more massive stars at the expense of 



the less massive ones (Bonnell et al. 1997 19981. However, by now we know that 



the presence of lower-mass companions can significantly reduce the growth rate of 



massive stars (Peters et al. 2010b). This is because neighboring stars will swallow 



material that would otherwise end up being accreted by the massive stars in the 
center. We call this process fragmentation-induced starvation and discuss it in 
more detail in Section |6.3| below. In reality both effects are present and influence 
the resulting stellar mass spectrum. 

Once a star has reached a mass of ~ 10 M Q it begins to ionize its environment. 
It carves out a bubble of hot and tenuous gas, which eventually will expand and 
enclose the entire stellar cluster. At this point no new stars can form and stellar 
birth has come to an end. We can observe the young cluster at infrared or even 
optical wavelength, as illustrated in Figure [6] 



6.2 Models of the Origin of the IMF 

There are three dominant schools of thought regarding the origin on the IMF each 
focusing on different aspects of gravitational collapse in the turbulent ISM. We 
point out, that the boundaries between these pictures are not clearly defined and 
that numerous hybrid models have been proposed. 

Core accretion. This model takes as its starting point the striking similar- 
ity between the shape of the observed core mass distribution and the IMF. This 
model assumes a one-to-one relation between the two distributions, such that the 
observed cores cores are the direct progenitors of individual stars or binary sys- 
tems. The factor of ^3 decrease in mass between cores and stars is thought to 
be the result of feedback processes, mostly protostellar outflows, that eject a fixed 
fraction of the mass in a core rather than letting it accrete onto the star ( |Matznerj 
and McKee 2000 ) . This model reduces the problem of the origin of the IMF to 



understanding the mass spectrum of bound cores. Arguments to explain the core 



mass distribution generally rely on the statistical properties of turbulence (Klessen 



|2001a||Padoan and Nordlund[ |2002l |Padoan et aT| [2007] [Hennebelle and Chabrier 
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2008 2009 



Chabrier and Hennebelle 


2010 


Veltchev et al. 


2011 



structures with a pure powerlaw mass spectrum. The thermal Jeans mass in the 
cloud then imposes the flattening and turn-down in the observed mass spectrum. 

Competitive accretion. A second line reasoning accounts for the fact that stars 
almost always form in clusters, where the interaction between protostars, and 
between a protostellar population and the gas cloud around it may become im- 
portant (Bonnell et al. 2001a|b Bonnell and Bate 



and Bonnell 2005 Bonnell et al. 2006 2008). In this picture the origin of the 



2002 


Bate et al. 


2003 


Bate 



peak in the IMF is much the same as in the core accretion model, it is set by the 
Jeans mass in the prestellar gas cloud. However, rather than fragmentation in the 
gas phase producing a spectrum of core masses, each of which collapses down to 
a single star or star system, in the competitive accretion model all gas fragments 
down to roughly the Jeans mass. Prompt fragmentation therefore creates a mass 
function that lacks the powerlaw tail at high masses that we observe in the stellar 
mass function. This part of the distribution forms via a second phase in which 
Jeans mass-protostars compete for gas in the center of a dense cluster. The cluster 
potential channels mass toward the center, so stars that remain in the center grow 
to large masses, while those that are ejected from the cluster center by TV-body 



interactions remain low mass ( Klessen and Burkcrt 2000 Bonnell et al. 2004 1 . In 



this model, the apparent similarity between the core and stellar mass functions is 
an illusion, because the observed cores do not correspond to gravitationally bound 



structures that will collapse to stars (Clark and Bonnell 2006 Smith et al. 2008) 



Importance of the thermodynamic behavior of the gas. One potential drawback 
to both the core accretion and competitive accretion models is that they rely 
on the Jeans mass to determine the peak of the IMF, but leave unanswered the 
question of how to compute it. This question is subtle because molecular clouds 
are nearly isothermal, but they contain a very wide range of densities, and it is 
unclear which density should be used. A promising idea to resolve this question, 
which is the basis for a third model of the IMF, focuses on the thermodynamic 
properties of the gas. The amount of fragmentation occurring during gravitational 



collapse depends on the compressibility of the gas ( Li et al. 2003 1 . For polytropic 



indices 7 < 1, turbulent compressions cause large density enhancements in which 
the Jeans mass falls substantially, allowing many fragments to collapse. Only a few 
massive fragments get compressed strongly enough to collapse in less compressible 
gas though. In real molecular gas, the compressibility varies as the opacity and 
radiative heating increase. Larson (2005) noted that the thermal coupling of the 
gas to the dust at densities n cr it ~ 10 5 cm -3 leads to a shift from an adiabatic index 
of 7 ~ 0.7 to 1.1 as the density increases above n cr u. The Jeans mass evaluated at 
the temperature and density where this shift occurs sets a mass scale for the peak 
of the IMF. The apparent universality of the IMF in the Milky Way and nearby 
galaxies may be caused by the insensitivity of the dust temperature on the intensity 
of the interstellar radiation field (Elmegreen et al. 2008). Not only does this 



mechanism set the peak mass, but also appears to produce a power-law distribution 



of masses at the high-mass end comparable to the observed distribution ( Jappsen 
eFaLl [20051). 
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Caveats. Each of these models has potential problems. In the core accretion 
picture, hydrodynamic simulations seem to indicate that massive cores should 



Clark and Bonnell 


2006 


Bonnell and Bate 


2006) 



The hydrodynamic simulations 
often suffer from over-fragmentation because they do not include radiative feedback 



from embedded massive stars (Krumholz 2006 Krumholz et al. 2007a Krumholz 



and McKee 2008 ) . Modeling the formation of massive stars consistently is a very 



active and challenging field of modern astrophysical research and discussed in more 
detail in Section |6.3| below. The suggestion of a one-to-one mapping between the 
observed clumps and the final IMF is subject to strong debate, too. Many of 
the prestellar cores discussed in Se ction [i] appear to be stable entities ( Johnstone 



et al. 2000 2001 2006 Lada et al. 2008 1, and thus are unlikely to be in a state of 



active star formation. In addition, the simple interpretation that one core forms 
on average one star, and that all cores contain the same number of thermal Jeans 



masses, leads to a timescale problem (Clark et al. 2007) that requires differences 



in the core mass function and the IMF. The criticism regarding neglect of radiative 
feedback effects also applies to the gas thermodynamic idea. The assumed cooling 
curves typically ignore the influence of protostellar radiation on the temperature 



of the gas. which simulations show can reduce fragmentation (Krumholz et al. 



2007a I . The competitive accretion picture has also been challenged, on the grounds 



that the kinematic structure observed in star-forming regions is inconsistent with 
the idea that protostars have time to interact with one another strongly before 
they completely accrete their parent cores (Andre et al. 2007). 



6.3 Massive Star Formation 



Because their formation time is short, of order of 10 5 yr, and because they build 
up deeply embedded in massive cloud cores, very little is known about the initial 
and environmental conditions of high-mass stellar birth. In general high-mass 
star forming regions are characterized by more extreme physical conditions than 
their low-mass counterparts, containing cores of size, mass, and velocity dispersion 
roughly an order of magnitude larger than those of cores in low-mass regions (e.g. 
Jijina et al.] p99| [Uaray and Lizano] [19991 |Kurtz et aL| [20001 |Beuther et U] 



2007 



Motte et al. 20081. Typical sizes of cluster-forming clumps are ~ lpc, 



they have mean densities of n ~ 10° cm , masses of ~ 10 M Q and above, and 
velocity dispersions ranging between 1.5 and 4kms . Whenever observed with 
high resolution, these clumps break up in even denser cores, that are believed to 
be the immediate precursors of single or gravitationally bound multiple massive 
protostars. 



Massive stars usually form as members of multiple stellar systems (Ho and 



Haschick 




1981 




Lada 


2006 


Zinnecker and Yorke 




2007 


) which themselves are 


parts of larger 


clusters ( 


Lada and Lada 


2003 


de Wit et al. 2004 


Testi et al. 



1997). This fact adds additional challenges to the interpretation of observational 



data from high-mass star forming regions as it is difficult to disentangle mutual 



dynamical interactions from the influence of individual stars (e.g. Goto et al. 
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2006 Linz et al. 2005). Furthermore, high-mass stars reach the main sequence 



while still accreting. Their Kelvin-Hclmholtz pre-main sequence contraction time 
is considerably shorter than their accretion time. Once a star has reached a mass 
of about 10 M Q its spectrum becomes UV dominated and it begins to ionize its 
environment. This means that accretion as well as ionizing and non- ionizing ra- 



diation needs to be considered in concert (Keto 2002 Keto 2003 2007 Peters 



et al. 2010a I. It has been realized decades ago that in simple 1-dimensional 



collapse models the outward radiation force on the accreting material should be 



significantly stronger than the inward pull of gravity (Larson and Starrfield 1971 



Kahn 1974 Wolfire and Cassinelli 1987) in particular when taking dust opacities 
into account. Since there are observations of stars up to ~ 100 M (Bonanos et al. 



2004 Figer 2005 Rauw et al. 2005 1 a simple spherically symmetric approach to 



high-mass star formation must fail. 

Consequently, two different models for massive star formation have been pro- 
posed. The first one takes advantage of the fact that high-mass stars always form 
as members of stellar clusters. If the central density in the cluster is high enough, 
there is a chance that low-mass protostars collide and so successively build up 
more massive objects (Bonnell et al. 1998). As the radii of protostars usually 



are considerably larger than the radii of main sequence stars in the same mass 
range this could be a viable option. However, the stellar densities required to 



produce massive stars by collisions are extremely high (Baumgardt and Klessen 



2011 ) and seem inconsistent with the observed values of Galactic star clusters (e.g. 



Portegies Zwart et al. 2010 and references therein). 



An alternative approach is to argue that high-mass stars form like low-mass 
stars by accretion of ambient gas that goes through a rotationally supported disk 
caused by angular momentum conservation. Indeed such disk structures are ob- 
served around a number of high-mass protostars ( Chini et al. 20041 12006 , Jiang 



et al. 


2008 


Davies et al. 


2010) 



that might have been present in the initial cloud and thus solves the opacity prob- 
lem. Radiation tends to escape along the polar axis, while matter is transported 
inwards through parts of the equitorial plane shielded by the disk. Hydrodynamic 
simulations in two and three dimensions focusing on the transport of non-ionizing 



radiation strongly support this picture (Yorke and Sonnhalter 2002 Krumholz 



et al. 2009 Kuiper et al. 2010 2011 ). The same holds when taking the effects of 



ionizing radiation into account ( Peters et al. 2010a|c|b 2011). If the disk becomes 
gravitationally unstable, material flows along dense, opaque filaments whereas the 
radiation escapes through optically thin channels in and above the disk. Even ion- 
ized material can be accreted, if the accretion flow is strong enough. H n regions 
are gravitationally trapped at that stage, but soon begin to rapidly fluctuate be- 



tween trapped and extended states, in agreement with observations ( Peters et al. 



2010a Galvan-Madrid et al. 2011 ). Over time, the same ultracompact Hn region 



can expand anisotropically, contract again, and take on any of the observed mor- 



phological classes (Wood and Churchwell 1989 Kurtz et al. 1994 Peters et al. 



2010c). In their extended phases, expanding Hn regions drive bipolar neutral 



outflows characteristic of high-mass star formation (Peters et al. 2010a) 
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Another key fact that any theory of massive star formation must account for 
is the apparent presence of an upper mass limit at around 100— 150 M Q (Massey 
2003). It holds for the Galactic field, however, it is also true for young star 



clusters that are massive enough so that purely random sampling of the initial 
mass function (IMF) (Kroupa 2002 Chabrier 2003) without upper mass limit 



should have yielded stars above 150 M Q (Weidner and Kroupa 2004 Figer| 2005 
Oey and Clarke 2 005[ | Weidner et al. 2010[ see however, Selman and Melnick 



2008). This immediately raises the question of what is the physical origin of 
this apparent mass limit. It has been speculated before that radiative stellar 
feedback might be responsible for this limit (for a detailed discussion see, e.g., 



Zinnecker and Yorke 2007) or alternatively that the internal stability limit of 



stars with non-zero metallicity lies in this mass regime ( Appenzeller 1970b|a 1987 



Baraffe et al. 2001). However, fragmentation could also limit protostellar mass 
growth, as suggested by the numerical simulations of Peters et al. (2010b). The 



likelihood of fragmentation to occur and the number of fragments to form depends 
sensitively on the physical conditions in the star-forming cloud and its initial and 
environmental parameters (see, e.g., Girichidis et al. 20111. Understanding the 



build-up of massive stars, therefore, requires detailed knowledge about the physical 
processes that initiate and regulate the formation and dynamical evolution of the 



molecular clouds these stars form in ( Vazquez-Semadeni et al. 2009) 



Peters et al. (2010b 20111 argue that ionizing radiation, just like its non- 



ionizing, lower-energy counterpart, cannot shut off the accretion flow onto massive 
stars. Instead it is the dynamical processes in the gravitationally unstable ac- 
cretion flow that inevitably occurs during the collapse of high-mass cloud cores 
that control the mass growth of individual protostars. Accretion onto the cen- 
tral star is shut off by the fragmentation of the disk and the formation of lower- 
mass companions which intercept inward moving material. They call this process 
fragmentation-induced starvation and show that it occurs unavoidably in regions 
of high-mass star formation where the mass flow onto the disk exceeds the inward 
transport of matter due to viscosity only and thus renders the disk unstable to 
fragmentation. 

As a side note, it is interesting to speculate that fragmentation-induced starva- 
tion is important not only for present-day star formation but also in the primordial 
universe during the formation of metal- free Population III stars. Consequently, we 
expect these stars to be in binary or small number multiple systems and to be of 



lower mass than usually inferred (Abel et al. 2002 Bromm et al. 2009). Indeed, 



current numerical simulations provide the first hints that this might be the case 



(Clark et al. 2011 Greif et al. 2011). 
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